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THE BLOW-UP OF ELECTROMAGNETIC FIELDS IN 
3-DIMENSIONAL INVISIBILITY CLOAKING FOR MAXWELL’S 

EQUATIONS 

MATTI LASSAS AND TING ZHOU 


Abstract: Transformation optics constructions have allowed the design of cloak¬ 
ing devices that steer electromagnetic, acoustic and quantum waves around a region 
without penetrating it, so that this region is hidden from external observations. The 
proposed material parameters are anisotropic, and singular at the interface between 
the cloaked region and the cloaking device. The presence of these singularities causes 
various mathematical problems and physical effects on the interface surface. In this 
paper, we analyze the 3-dimensional cloaking for Maxwell’s equations when there 
are sources or sinks present inside the cloaked region. In particular, we consider 
nonsingular approximate invisibility cloaks based on the truncation of the singular 
transformations. We analyze the limit of solutions when the approximate cloaking 
approaches the ideal cloaking in the sense of distributions. We show that the so¬ 
lutions in the approximate cloaks converge to a distribution that contains Dirac’s 
delta distribution supported on the interface surface. In particular, this implies that 
the limit of solutions are not measurable functions, making them outside of those 
classes of functions that have earlier been used in the models of the ideal invisibil¬ 
ity cloaks. Also, we give a rigorous meaning for the “extraordinary surface voltage 
effect” considered in physical literature of invisibility cloaks. 

Key words: Invisibility cloaking, Maxwell’s equations, transformation optics. 

1. Introduction 

The transformation optics-based cloaking design, since proposed in 2006 in [19, 
31], has attracted most attention among many theoretical proposals for realizing 
invisibility, with widely reported experiments and data accumulated due to the de¬ 
velopment of metamaterials. The basic idea is that, the coordinate transformation 
invariance for certain systems, such as those describing electromagnetic/optic or 
acoustic wave propagations, makes modifying the background medium of a local 
region Id in a certain way undetectable when observe far away. In particular, when 
singular spatial transformations are used, like the one that blows up a point into 
a small bounded domain D, this undetectable customization of medium (known as 
the push-forward by the singular transformation) can have the structure of a fixed 
layer I1\D (the cloaking device) surrounding “arbitrary” media in D (the cloaked 
region). Namely, whatever in D along with the cloaking device is invisible! Such 
singular coordinate transformations to create a “hidden pocket” surrounded by a 
layer of degenerate material, were proposed in 2003 in [14] for electrostatics as 
a counter-example of uniqueness of an inverse problem, known as the Calderon’s 
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problem. In [19, 31] a customized cloaking layer was proposed to be created of meta¬ 
materials. Such layer should be designed so that it bends the light rays or detour 
the electromagnetic waves away from the cloaked region and return them back, as if 
they penetrate the region straightly, making observation outside indistinguishable 
from that of the empty background space. 

The main difficulty in analyzing this scheme rigorously lies on the prescribed sin¬ 
gular (non-regular) medium in the cloaking layer, see discussion in [3, 9, 10]. To be 
more precise, we consider in this paper the electromagnetic cloaking, where the pre¬ 
scribed medium in the cloaking layer Q\D has electric permittivity e(a:), magnetic 
permeability Since they are the push-forward of the ambient medium param¬ 

eters by a singular transformation (detailed in Section 2), one of the eigenvalues 
of e{x) and ]l{x) degenerates to zero at the exterior of the cloaking interface dD'^. 
Suppose D is filled with an arbitrary medium {gii{x)^Si{x)) to be cloaked. Let us 
consider the waves propagating in domain C consisting of medium whose per¬ 
mittivity and permeability are given by (^, e) in n\D (the cloaking layer/device) 
and (/ii,£i) in D (the cloaked region). Such cloaking gives counterexamples for 
uniqueness of the inverse problems for Maxwell’s equations. As we want to consider 
also cloaking of active bodies, we add a current source J in the domain D. To study 
this, we will consider the case when /ri(x) = /tq and £i(a;) = Eq are constants. We 
note that the results of this paper can with small modifications be extended to the 
case where fJ-i(x) — p,o and £i(x) — eg are compactly supported functions in £>. 

Since (/i, £) are singular near dD, the key mathematical question is, in what 
sense the solutions to Maxwell’s equations 


curl E — iujB = 0, curlll-h icuD = J, on fl\DUD, 
exist. Here the constitutive relations are given by 




on fl\D, 
on D. 


on fl\D, 
on D. 


Moreover, one would like to specify what kind of boundary conditions should ap¬ 
pear at the interface dD, where the singularity happens. There have been many 
different proposals to model the above Maxwell’s equations with singular coeffi¬ 
cients and different boundary conditions have been proposed, depending on what 
spaces the solutions of the equations are assumed to belong. In the present paper 
our aim is to consider approximative cloaks and consider what is the limit of the 
solutions as it approaches an ideal cloaking. Before formulating our results, let us 
review the earlier proposed models. 

In [7], the concept of finite energy solution (FES) was formulated. These are 
distributional solutions to Maxwell’s equations for which all physical fields, that is, 
the electric and magnetic fields E and H as well as the electromagnetic fluxes D and 
B are C^-valued measurable functions on H. The “tilde” here refers to the fact that 
these solutions satisfy Maxwell’s equation with the cloaking material parameters £ 
and Jl. Moreover, the energy norm with degenerate weight is required to be finite, 
that is 
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where, as everywhere below, we have used the Einstein summating convention by 
summing over indecies j and k that appear both as sub- and superindecies. It is 
shown in [7] that a hidden boundary condition must be satisfied by the FES at the 
interior of the interface, namely 

X E\dD- =0) X H\gD- = 0. 

These, physically known as PEC and PMC conditions, form an over-determined set 
of boundary conditions for Maxwell’s equations in D, and could be only possibly 
satisfied for cloaking passive media (i.e., J = 0 or J is a non-radiating source that 
would produce in the free space a compactly supported electromagnetic field). 
But for generic J, including ones arbitrarily close to 0, there is no FES-solution. 

In [33, 34], solutions with various boundary conditions at dD^ are considered 
in the context of self-adjoint extensions of the Maxwell’s operator, that are com¬ 
patible with energy conservation. There, the fields E and H are assumed to be 
measurable functions that are in the domain of the closed quadratic form, as well 
as in the domain of the corresponding self-adjoint operator. We call these the op¬ 
erator theoretic (OT) solutions. It is shown there that the Maxwell operator A 
with degenerate coefficients s(x) and /l(x) is essentially self-adjoint and its unique 
self-adjoint extension has a domain V(A) C such that {E,H) G 'D(A) have 

one-sided traces (from the outside of D) satisfying 

V X E\Qjg+ =0, V X El\Qjg+ = 0, 

due to the degeneracy of parameters at dD'^. With regular medium to be cloaked 
in D, functions (E, iJ) G 'D(A) have also one-side traces, from the inside of H, 
satisfying 

(1) v ■ (curl H)\qd- = V ■ (curl = 0. 

Notice that this is from the point of view that self-adjoint extensions in the domain 
n are direct sums of those in ^A\D and D, hence the solutions are considered com¬ 
pletely decoupled from each other. 

In [8, 12, 13, 15, 18, 21, 25, 26, 27, 28], for the purpose of both physical realization 
and analysis of the waves behavior, the regularization of the singular scheme is 
studied as the fundamentals of approximate cloaking design. In this paper, we 
will consider the same approximation scheme as in [21] for electromagnetic waves. 
The detailed non-singular approximation will be introduced in Section 2. Roughly 
speaking, instead of using the transformation that blows up a point into the cloaked 
region D — Bi where Br C is the ball of radius R centred at the origin, in this 
idealized spherical geometric setting, we use a non-singular transformation that 
blows up a small ball Bp with radius p (the regularization parameter) into Bi. 
The transformation medium is then regular for fixed p > 0 and appears as a small 
inhomogeneity in the empty space. The well-understood solutions for regular media 
provide a tool to scrutinize the extreme case, that is taking the limit p —>■ 0. 

As seen above, there are different, partly contradicting, alternatives for rigorous 
models of ideal invisibility cloaks. As physical attempts to build up an invisibility 
cloak from metamaterial are always based on approximate constructions, we pro¬ 
pose here a change of the point of view for defining a model of an ideal invisibility 
cloak: We consider a limit of the solutions in an approximate cloak as p —>■ 0, 
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in the sense ofdistributions (i.e., generalized functions), and describe the plausi¬ 
ble solutions in an ideal cloak as these limits of the solutions. We show that the 
solutions in the approximative cloaking structures converge to a distribution that 
contains Dirac’s delta distribution supported on the interface dD. This implies that 
the set of measurable functions might be too small as the domain of the singular 
Maxwell’s operator, when considered as the limit of electromagnetic approximate 
cloaking. We note that the ideal and the approximate cloaks we study here have 
no energy absorption, that is, the conductivity is zero. 

The results of the paper can be summarised as follows: Let Ep and Hp denote 
the electric and the magnetic fields in an approximative cloak in D as described in 
(15), satisfying the boundary condition v x Ep\dB 2 = /■ We consider the limit, in 
the sense of distributions, En = lim Eg and Hn = lim El^ of the solutions of the 

p-fO p-fO 

Maxwell equations, that is, the limit when the approximative cloak approaches a 
perfect cloak. Then 

(1) We show that the limit {Eq, Hq) is not in the class of finite energy solutions. 
This is compatible with the non-existence result of finite energy solutions 
given in [7]. 

(2) We show that the limit {Eq,Hq) is not a measureble function and thus 
it does not belong in the class of operator theoretic solutions studied in 
[33]. However, if {Eq,Hq) is decomposed to a sum of a measurable func¬ 
tion [E^^Hm) and a delta-distribution supported on the surface dD, the 
function part {Em,Hm) satisfies the non-standard boundary conditions (1) 
on dD described in [33], that is, certain traces of {Em,Hm) are zero. Sur¬ 
prisingly, we observe that the corresponding traces of the fields Ep and Hp 
on dD do not converge to zero, see Remark 4.6. 

(3) The convergence of the solutions in the so-called virtual space and physical 
space are different: In the physical space the limit of {Ep,Hp) contains a 
delta-distribution part but in the virtual space, images of the helds given 
in (17), converge in the L^-sense as described in [21]. 

(4) The blow up of the wave, that is, the appearance of the delta-distribution 
in the limit p —>■ 0 happens with all frequencies. Thus the observed blow 
up is different to the destruction of the cloaking effect appearing at the 
eigenfrequencies of the inside of the cloak studied in [8, 12, 15]. 

(5) The high concentration of waves on dD has been observed in the physical 

literature by Zhang et. al. [36], and is called the “extraordinary surface 
voltage effect”. We give a rigorous formulation, as a convergence result in 
the weak topology of of this effect. 

We note that in this paper we consider only the most commonly used approxima¬ 
tive cloaks that are equivalent to multiplying the singular electric permittivity e{x) 
and magnetic permeability ^(x) by a characteristic function. Such approximation 
could be implemented also in different ways, see [2, 4, 5, 8, 15, 18, 17, 16, 21, 25, 
26, 27, 28], and an interesting question is how the delta-distribution part of the 
limit solution {Eq,Hq) depends on how the approximation is used. However, this 
question is outside the scope of this paper. 

The rest of this paper is organized as following. In Section 2, we describe in 
details the transformation optics-based ideal cloaking and the regularized approx¬ 
imate cloaking in the spherical geometric setting. The main result is formulated 
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in Theorem 2.1 and the proof is presented in Section 4, where the Dirac’s delta 
distribution is proved rigorously to appear at the exterior of the interface dD^, as 
a limit of the high concentrated electromagnetic waves in the approximate cloak. 
In Section 3, an example is presented to demonstrate the physical motivation and 
insight of the existence of such condition. Our estimates rely heavily on the spheri¬ 
cal harmonics and asymptotic of Bessel and Hankel functions, of which we provide 
some basics in the Appendix. 

2. Singular ideal electromagnetic cloaking and the regularization 

How the transformation invariance of the system amounts to the undetectability 
of certain non-ambient media with respect to the observation is best seen when 
formulated as inverse boundary value problems for underlying PDEs (for near- 
held measurements) or inverse scattering problems (for far-held measurements). 
This is also why the layer-structured medium for cloaking was hrst studied in 
[14] for electrostatics as a counter-example of uniqueness of an inverse problem, 
known as the Calderon’s problem. Here we consider a similar construction for 
electromagnetism, that is for Maxwell’s equations. In particular, we prescribe time- 
harmonic incident waves with time frequency w, which results in the time harmonic 
equations as following. 

Let H denote a bounded domain in with smooth boundary. Consider Maxwell’s 
equations for time-harmonic electric and magnetic helds (A, H), viewed as 1-forms, 

(2) V X A = iojB, V X iJ = —iujD + J in H 

where the 2-forms B, D and J denote the magnetic induction, the electric displace¬ 
ment and the current density respectively. Furthermore, we specify the constitutive 
relations as 

(3) D = eE, B = nH, 

where e and /i represent the permittivity and permeability of the material in H . 

In the following we assume that e, ^ are in L°°(S1)^^^, and satisfy 

(4) < C^£:(ai)C < CmICP, < ^^^^(a:)^ < CmICP 

for some constants Cm,CM > 0 and all a: S H and ^ G K^\{0}, We remark that 

(4) are physical conditions for regular EM media. Given J G Lp'iyiY’, we denote by 
a subset of (Div; dU) x (Div; dO) given by 

'■= |(u X E\qq,v X i?|an) \ E,E[ G iL(curl; H) satisfy (2) and (3) with j| 

where 

7L(curl; H) := {u G L^(H)^; curl u G L^(H)^}, 

and 

(Div; 3H) := {f G E[~^ (9H)^; i ■ v = 0 a.e. on and Divf G 7L“^(9H)} 

with Div denoting the surface divergence on 9H. Also, we denote := j 
with J = 0. 

The set j is known as the Cauchy data set which encodes the full exterior 
(boundary) measurements of electromagnetic fields. The inverse problem is then to 
understand the dependence of j on the parameters in order to recover 

the latter from the former. Knowing this, the design of invisibility is to find the 
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parameters to be prescribed in a cloaking device such that its Cauchy data is in¬ 
distinguishable from that of the vacuum background, independent of the object to 
be cloaked in the device. 

First, let us define M~'^ = for an matrix M. Consider a transforma¬ 

tion X = F{y) : n —n between two bounded domains fl, n C with smooth 
boundaries. Assume F is bi-Lipschitz and orientation-preserving, and denote by 
M := DF{y) = the Jacobian matrix of F. The pull-back fields of the 

solution {E,H) G iJ(curl;n) x iJ(curl;fJ) to (2) by F~^, are given in 17 by 

Fix) ={F-^rE{x) := {M-^E) o F-\x), 

(5) _ 

H{x) ={E-^)*H{x) := o F-\x), 

(6) J{x) = (F-i)* J(x) = ([det(M)]-^MJ) o E-\x). 

Then we have that {E,H) G iJ(curl;17) x iJ(curl;17) satisfies Maxwell’s equations 

(7) V X E = iojJlFl, V X = —iojeE + J in 17, 

where V x denotes the curl in the x-coordinates, and e, J1 are the push-forwards of 
e, /i via F, defined by 

(8) e{x) = F^e{x) := ([det(M)]“^M • e • M^) o F~^{x), 

and similarly for 'll = F^y,. Moreover, if we assume Ajan = Id, using Green’s 
identity, it is directly verified that 

nu) _ nu) 

This summarizes the basics of transformation optics in a rather general setting, 
which we shall make essential use of in the construction of singular ideal cloaking 
and regularized approximate cloaking. Throughout the paper, we will focus on a 
spherical geometry of design. 

Next, we consider singular coordinate transformations. Consider the map 

(9) Fi : B2\{0} ^ iJ2\F, F’i(y) = 0 < |y| < 2 

which blows up the origin to Bi while keeping the boundary dB 2 fixed. In the 
cloaking layer x G B 2 \Bi, we prescribe the EM material parameters given by 

(10) yix) = £(x) = {Ful){x) := 

|det(77Fi)| 

where I is the identity matrix, representing the homogeneous vacuum background 
material. In the region Bi to be cloaked, we consider an arbitrary but regular EM 
medium (eoiMo) satisfying (4), i.e., '^{x) = yo{x) and e'(x) = eo(x) for x G Bi, 
which can be viewed as the push-forwards of {yo, £o) in Fi by F 2 = Id. We denote 
the “glued” transformation by 

(11) F=(Fi,F2): (B2\{0},Bi)^(B2\B;,Bi). 

Noticing that Fi is a radial dilation and by some simple calculations, we have that 

(12) /I(x) = £(x) = 2 ^^^j e^-|-2ee, 1 < |x| < 2, 

X 
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where and are respectively, the unit projections along radial and angular 
directions, i.e., e^. = / — xSF^ eg = x = xl\x\. It is readily seen that as one 
approaches the cloaking interface the medium in the cloaking device becomes 
singular, in the sense that e and Ji no longer satisfy the condition (4) (the eigenvalue 
along the radial direction degenerates). 


2.1. Construction of regularized approximate cloaking. For approximate 
acoustic cloaking by regularization, Kohn et ah, in [15], proposed blowing up a 
small ball Bp to Bi using a nonsingular transformation Fp which degenerates to 
the singular transformation Fi in (9) as p —>■ 0, while Greenleaf et ah, in [8], 
proposed truncating the singular medium in (12) to for i? > 1. For the 

present study, we shall focus on the ‘blow-up-i3p-to-i?i’ regularization. 

Let 0 < p < 1 denote a regularizing parameter and set 


(13) 


a = 


2(1-P) 

2-p 


b = 


1 


P' 


Consider the nonsingular transformation from B 2 to B 2 defined by 

= (a + &|2/|)|fy, for p < jpj < 2, 
for ||y| < p. 


(14) 


a; := Fp{y) = 


F^^\y) = % 


Our approximate cloaking device is obtained by the push-forward of a homogeneous 
medium in B 2 \Bp by Fp^\ Suppose we hide a regular but arbitrary uniform EM 
medium (eo,po) in the cloaked region Bi. Then the corresponding EM material 
parameter in B 2 in the physical space is 


(15) 


(ep{x),Hp{x)) = 


(^Oj Po)i 


for 1 < |a:| < 2, 
for |a:| < 1. 


The EM fields {Ep,Hp) G H{cml-B 2 ) x iL(curl; 52 ) corresponding to {B 2 ;ep,fip} 
then satisfy Maxwell’s equations 


(16) 


W X Ep = iujfipHp, V X Hp = —iuSpEp + J in i? 2 , 
V X EpldB:, = / G iL“fo^(Div;5B2) 


where J G (L^(Bi))^ is the current source. At this point, we assume that it is 
supported in a smaller ball Br^ with radius 0 < ri < 1. Then the pull-back EM 
fields 


(17) (Ep, Hp) = {{Fp)*Ep, {Ep)*Hp) G iL(curl; B 2 ) x i7(curl; B 2 ) 


satisfy Maxwell’s equations in the virtual space of y with parameters 


(ep(2/),Pp(2/)) 


{I.I) P<|2/I<2, 

((Ei"))*£o,(i"i'Vpo) |y|<P 


and 

J=(Ff )V. 

Moreover, the observation CG of the whole cloaking object {B 2 ; ?p, Pp} is shown 
to be identical to , that of a small inhomogeneity of radius p in the background 
vacuum space. In particular, the regularized cloaking is expected to converge to the 
ideal cloaking as p shrinks. The order of such convergence was discussed in [15, 21] 
for both Helmholtz equations and time-harmonic Maxwell’s equations. Here we 
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reproduce some calculations from [21] for Maxwell’s equations in order to pursue 
our discussion. 

Due to the two-layered structure, the following notations for EM fields will be 
adopted 


~ _ f E+{x), iovx€B2\Bi, 

^ ' I E-{x), for xe Bi, 

in the physical space and 

p f E+{y), toryeB2\Bp, 

^ ' 1 (y)> for y G Bp, 



H+{x), for X e B 2 \Bi, 
Hp{x), for a; G i?i, 


f H+{y), forj/GB 2 \Bp, 

1 for y G Bp, 


in the virtual space, and the fields satisfy the following transmission problems 


(18) 


and 


(19) 


{ V X E+ = iujy.p{x)H+, V2< H+ = -*wep(a;)E+ in B 2 \Bi, 
V X Ep = iu}fj,oE[p , V X Efp = —iuje^Ep + J \n Bi, 

^ ^ ^\dB+ = ^ ^ ^p loB”’ ^ ^ ^p'l9B+ = ^ ^ ^p laB“’ 

1/ X E+\e,B2 = /■ 

{ V X £1+ = iujH+, V X £f+ = -iujE+ in B 2 \Bp, 

W X Ep = iuJHp{y)H-, V x i?" = -iuj£p{y)E- +J in Bp, 
^ ^ '^P~laB+ = ^ ^ ^p Usy’ ^ ^ -^p'Ub+ = ^ ^ Hp IdB-’ 

V X E+|aB2 = /■ 


Now we are ready to present our main theorem. 


Theorem 2.1. Let E and E[ be 1-forms satisfying the Maxwell’s equations on 
^2\{0}, at frequency uj > 0 which is not an eigenvalue for the background Maxwell’s 
operator, 


( 20 ) 


r V X E = iu>E[, V X i£ = —iojE on B 2 

\ B X E\qb^ = f, 


and let and Hq be solutions to 

J V X Eq = zw/iQi^Q, V X Hq = — iujEoEQ J 
\ V ■ EqIobi = V ■ £lo|aBi = 0 


( 21 ) 


on Bi 


where J is supported on B^ for some ri < 1. Moreover, suppose that E and El 
belong to L^{B 2 ',M.^) such that 


( 22 ) 


{E,H) 


{Ei^E,Ei^H), in B 2 \Bi 
{E o,Ho), in Bi 


where Ei is the singular transformation given by (9). Then we have 

(23) Ep^E + a[J]6aB, Hp ^ H + I3[J]5 ob^ as p ^ D 

in the weak topology of E[~^{B 2 ), where 5aBi is the Dirac’s delta distribution sup¬ 
ported on dBi and a[J], (3[J] are smooth and depending on the source term J (see 
Remark 4-5)- Moreover, we have (see Remark 4.2) 


(24) 


{curl £f)| 9 B- = B • {curl E)]^^- 


= 0 


and 

(25) 
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^ ^ ^\dB+ — ^ ^ ^\dB+ ~ 0 - 


We note that despite (25) 


(26) 


hi := lim v x Hp\oBi and := lim i' x -Eplasi- 


may be non-zero. (See Remark 4.6) 

3. Motivation: A physical example of scattering from the half-space. 

Before the analysis of the waves in an approximative cloak, let us us consider as 
a motivation on a simple example of scattering of a plane wave from homogeneous 
half-space that resembles the approximate cloak near the cloaking surface. To 
begin with, we decompose the space to half spaces U+ = {{x,y, z); x> 0} and 
U- = {(x, y, z)] X < 0}, and their interface S = {(x, y, z); x = 0}. We assume that 
the electromagnetic parameters in [/_ correspond to vacuum, that is, e_ = J and 
/i_ = I. In {7+ the electromagnetic parameters are given by constant matrices. 


^t) M+= diag (/rL/iy 


£+ = diag {e+,e: 


X ^ 


Below, ex = (1,0,0), Cy = (0,1, 0), and Cz = (0,0,1) are the unit coordinate vectors. 

To study the case that is close to the 3-dimensional approximate cloak studied 
in this paper, we let 



To illustrate, we consider an incident plane wave whose magnetic field is parallel 
to y-axis. It comes from the domain U- to the interface S with non-zero incident 
angle and scatters and refracts at the interface. More precisely, we let magnetic 
field be Hp{x, z) = h{x, z)ey, where 


/i(x, z) = h+{x, z) := for x > 0 


and 


h{x, z) = h-(x, z) := -b ^">)ey for x < 0. 


Recall that this corresponds to the decomposition of the wave in U- into a sum 
of incident and reflected plane waves. At frequency w > 0, given the incident 
amplitude h™ S C and the ^-component of the incident wave vector k'^ G (0,a;), 
we consider the behavior of the transmitted wave in U- with different values of 

p > 0. 

The field Hp{x, z) satisfies Maxwell’s equations in U± if 



(27) 


This implies 



(28) 


Correspondingly, the electric field is of the form Ep = ExCx + EzCz where 
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Now the transmission conditions on S, 


fix X i?p|s-Ca; X 

Cx X -£^p|e“ ~ X -£/p|x;+ 

give /cj = and 

(29) = 

(30) = 

Since k~ = kf, equation (28) implies 

fc+ = ^et )2^ = ^^uj'^ - ■^{kt)'^ = ik+p~^{l + 0{p~^)). 

Note that when p —>■ 0, the expression under the square root becomes negative. We 
have choosen above the positive sign for the imaginary part of the square. Denote 
below t{p) = ImA:+. Then t{p) —)• oo as p —)■ 0. 

By solving (29) and (30) we obtain 


Thus 


h+ 





(l - 

\ ^kx kx J 


Hp{x,z) 

Hp{x,z) 


^-t(p)x+ik^ z- 


2kx + t{p)i 

k-x+ik-z _ M _ 


'VI 

AkZ 


2kx +t{p)i 


^—ik^ x-\-ik^ z \ ^ 


■y, 


for X < 0. 


Next, let Xb+(s) be the characteristic function of the set M+. Since 

lim XR+(s)^e“‘® = (5o(s) 
t—>-oo t 

in the sense of distributions in I?'(M), we see that 

lim Hp{x,z) = XK_(x)(e*'=^'"+*'=^'^ - - Aik-h^^6o{x)ey, 

p —>-0 

in the sense of distributions in I?'(IR^). This means that the magnetic field Hp 
tends to a generalized function that is a sum of a measurable function and a delta 
distribution as p —?► 0. Note that here the delta distribution component is caused 
by the concentration of the waves in a thin layer near the interface in the region [/_. 
Note that the blow up of the fields at the interface, that causes the delta distribution 
to appear, causes the boundary also to reflect the incoming wave perfectly with the 
reflection coefficient —1. Similar considerations to the above one, with slightly 
different setting where p is constant in the whole space, are done in [36]. Also, 
the physical interpretation of the blow-up of the fields at the interface as infinite 
polarization of the material is analyzed in [36] in detail. Our aim is to show that 
similar phenomenon appears in the 3-dimensional approximate cloak as the cloak 
tends to the ideal one and analyze the convergence of the electromagnetic fields in 
the sense of distributions. 
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4. Limiting behavior at the interface 


4.1. Spherical harmonic expansions. Through out the rest of the paper, we 
assume that the cloaked medium is uniform, i.e., Eq and /iq are constants. Then the 
EM fields and (i?+,iL+) both admit the spherical harmonic expansions, 

given by 

(31) 


E„ = e, 


- 1/2 


oo n 




n—1 m——n 

oo n 


E E + an V X 


n—1m——n 




for X e Bi\Bri, where k := {fioEoY^^, and 

OO n 

< = E E + CV X + CiV™. + CV X 


(32) 


n—1m——n 
oo n 


Hp =- E E + T^V X + CV X iV„- 

^ lUJ “ ’ , , , 


n—1 m——n 


for y G B 2 \Bp, where the basis vectors are defined in Appendix A by (67). Notice 
that the last two terms in the expansion of E~ with respect to andVxlV„-^ 

represent the radiating EM field generated by J. Hence, the coefficients p™ and 
are determined. 

In terms of the vector spherical harmonics (69), we have the expansion of the 
boundary condition 


OO n 


(33) 


/ = E E Snifi2u:r+fL^lv:r)- 


n—1 m——n 


where Sn '■= \/n{n + 1). By (74), (75) and the boundary condition v x E+ = fi 
we obtain 

(34) Chi^\2u;) + 7rjn(2w) = C ^n(2u;) + V^Jn{2co) = 2/^^. 

The transmission conditions in (18) reads 


(35) 

which gives 


X X E \gg+ — p{y X E la_B+) — x x E 


idB-y 


(36) 


PC^nH^P) + Plnjn{0jp) =Eq (tt™ (/cw) + (/cw)), 

ck^^Huiujp) + v^Mujp) (/3rJn(A:cc) + g™H„(fcca)). 


Similarly, the transmission condition on the magnetic field implies 

kc^Hniujp) + k'j'^Jniujp) =PQ^^^{a’^Jnikuj) + Pn'Hnikuj)), 

pdnh^nH^P) + PPnin{ujp) =Po (fc/3™ j„ (fcw) + kq^ hl^'> {ku})) . 


( 37 ) 
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Solving (36) and (37), we have 


(38) 


m I m , if _m m 

C-n — flln + tlPn : “ra 

4 . _m I if _m om 


I m , if rn 

hjn + f2Pn > 

I I if m 

tiVn + J 


where 


(39) 




^2 :=!=;— 


^3 


ti := 


1 

1 

Dn 

1 

1 


£Q^^'^kJn(ljjp)jn{kuj) - pjn{0Jp)Jn{kuj) , 

kpJn{Lup)h^^'> (ujp) - kpjn{ujp)'Hn{pJp) , 
PQ^^^kJ„{ujp)jn{kuj) - £Q^^‘^pjn{u}p)Jn{ku}) , 
pJniuJp)hi^\ujp) - pjn{ujp)Hn{ujp) ; 


and 


(40) 




to :=V^ 


:= — 


:= — 


1 

I^n 

1 

1 

1 


h^^\kuj)Jn(kuj)-'Hn{kuj)jn(koj) , 

kh^^'> {ku})'Hn{ujp) - p'Hn{ku})h^^\ujp) , 

Jn{kuj)h\l\kuj)-'Hn{kuj)3n{kuj) , 

PQ^^^kh^^\kuj)nniujp) - p'Hn{ku})h^^\ujp) , 


with 


Dn =Po ph^r!:H^P)Jn{kuj) - Eg ^ kHn{uj p)jn{kuj) , 
D'n phn\^P)Jnikuj) - Pg^^"^ kUniuJ p)j„ikuj) . 


Plugging into (34), we obtain 


(42) 


M-p'!;:t[h^P{2Lu) ^ _ 2f^± - t'^q^nn{2uj) 

tihn\2ljj) + jn{2uj) t3'Hn(2u;) + J'n{2u!) 


In this notes, we are interested in the effect of an active source cloaked in Bi, 
hence assuming no incident wave for simplicity, i.e., the case of zero boundary 
condition / = 0. Therefore, we have 


-p^t'^h^,^\2u2) ^ _ -t'^q^nn{2u2) 

tihn\2io) + jn{2ui) t3'Hn{2oj) + Jn{2uj) 
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We will also need the following asymptotic estimates derived from (77) for n> 0 
iTi{n -\- 1) fUJ\ 2n+l 


— 


t4 = 


(44) 


-(1 + Op.o(p)), 

r(?T, -I-3/2)/ig ' kuinjn{kui) 


T{n + l/2)r{n + 3/2)n 
{2n + l)v^ 


t3 = 


2 i^ Jn{kuj)hn\kui) —'Hn{ku))jn{kuj) 


(|)”^%”+Ml + Op-o(p)) 


r(n -I- 1/2)^q ^^^knjn{ku}) 

3n[kuj) 


where |Op_>.o(p)| < C(n)p for some constant C(n) depending on n. 

It suffices to show the proof for the electric field only, and that of magnetic field 
is obtained by symmetry. 

4.2. Normal components. Consider 

(45) j {x-Ep)(pdx = 

Jb2\b;^ 

for (fi G Hq{B 2 ), which admits the spherical expansion 

oo n 

(46) ip{x) = ‘Pn{\x\)Y;^{x), X e B2, 

and necessarily 


_ I _ {x ■ Ep )(pdx + / _ {x-E+)(pdx. 

B2\Sri JBi\Br-i JB2 \Bi 


n=0 m— — n 


oo n 


oo n 


\^n{r)\‘^rdr <oo, 


n—0 m——n '■ 


n—0 m——n * 


dr 


^nir) 


rdr < oo. 


In the physical region n < \x\ < 1, by (74) and (76), we have 

OO n ^ 

(48) x-E;= E E M \Pn3n{kix\x\) + h^^\ku:\x\)\ yr(i). 


n—1m——n ' ' 

Then the first integral of (45) is given by 


„ OO n 

/ _ {x-E-)ipdx= / E E ^™p(^ ^ 

3 BAB,., Jr, _-- 


IB,\B,^ 


’"I n—1 7n——r 


where 


[/3rjn(fcwf) 

Moreover, we can show 


Lemma 4.1. Suppose that J is supported in B^ for ri < I. Then we have 

P pi OO n 

I _ {x ■ Ep)(pdx = J E E 

’"I n—1 m— — r 


( 49 ) 


lim 

p—>0 


Bi\B,^ 
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The word “should” 
removed.- Matti 


Proof. Given that the current source J is supported on the radiation part 
(only depending on J) of E~ ^ that is, 


oo n 


(50) 




n—1m——n 


is C°° for I a; I > ri, hence is uniformly convergent for r G [ri, 1] and independent of 
p. Furthermore, we have 


(51) 


\Slqnh^n\ku}^r\ <M (1 + n) 


-M 


for r > ri 


for arbitrary M > 0, where A <m B represents A < CmB for a constant Cm 
only depending on M. Using (77) (or equivalently (78)), let iVi > 0 be such that 
Ni » kuj, then for n > Ni, \hn\kujr)\ is monotonically decreasing for r G [ri, 1], 
implying 

(52) 


\hn\kujri)\ 


By (38) and (43) we have 


am _ {ht'i - 

^n,p J. n / (rt. .\ I rr /o. 


^3^n(2w) + .7„(2a;) 


From (44), for n > 0 

— = -ITT- 


n + 1 


nr(n + l/2)F(n + 3/2) jn{kuj) 
which implies for n > 0 

am ^ {kuj) 


h^n\ku)ruj\^r,+ l 

TAi^\2) P (1 + Op^o(p)), 


jn{ku}) 


C(l + O,^o(p)). 


Therefore, 

(53) 


gm — lijn /3"" _ jkuj) 

^n,0 “ bn,p ~ „• q-n 


p—>-0 


jn{ku 


and the convergence is uniform. As a consequence, for p << 1 << 7Vi < n and 
r G [^ 1 ,1], the other term of can be bounded by 




< 


Pn,pjn{kuJ^ 


qnhn'’ {kwr) 

jn{kujrfhn\kuj) 




jn{kuj)hn\ku}r) 




<\Q^r 


'2-1/2 




Therefore, the series converges uniformly (independent of p) and by dominated 
convergence theorem the lemma is proved. □ 

Remark 4.2. The proof actually shows that 


lim / 

JBi\BA 


{x ■E)pdx= (x ■ Eq )p dx 
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where 


oo n 


X ■ (a;) := XI H [ l ^ n , ojn { kuj \ x \) + h^^\kuj\x\)\^ Y^{x). 

n—1m——n ' ' 

By the uniform convergence of (53) and estimate (52), we have 

* ■ ^0 \dB~ ~ 

Notice that the medium inside Bi is regular and J is supported away from dBi, 
this proves (24) in Theorem 2.1, which is consistent to Weder’s definition [33]. 


Now we focus on the second integral of (45). First of all, 

1 


(54) 


where 


/S2 \Bi 


{x ■ E'f)ip dx = 


B2\Bp b 


iy ■ K(Fpiy)) \DFp{y)\ dy 


p2 ^ ^ 


= j£Y. c.(e * 

P n—1 m——n 


{a + brY 


YYpY) ■■= Si dY:^phY\u:r) + T]Y,pjn{i^r) (pff{a + br) 

and := dff and := ly™ to indicate the p dependence. We separate the series 
into two parts: let iV 2 > 0 be such that max{2w, fcoj} << iVi < N 2 (as defined in 
(77)) and consider 

- /'E E dr. hip) := /' E E ^Zpir) dr. 

P n—1m——n 


’P n — N 2 -\- lFa——n 


First we have 

Lemma 4.3. 


As n 

(55) limA(p) = E E - 

n—1 m——n 


Jn{kuj)hn\kuj) - 'Hn{kuj)jn{kuj) 


Proof. We first consider the limit 

p2 -^2 n 


kn{n + l)jn{kui) 

(a + brY 




(56) 


lYffJ EE FlvY,pjn{uJr)ip^{a + br) 


dr 


P n—1 m——n 
N 2 n 


rv2 iL pZ 

= E E ^n'^YYo^Y,p jn{u}r)pff{a + br) 

- 1 ^ 'Jo 


{a + brY 


dr. 


n—1m——n 


From (43), (38) and (44), for p << 1 and 0 < n < iV 2 

-t^H„(2a;)C 


\Pn\ = 
Mnl = 


tsTdniYuj) + ffn(2uj) 

t'MMqY 


tsNni^Uj) + ffn(2uj) 


^N2 Idn IP 

|„m I „n+l 

N2 | 9 „ \P 


< 


( 57 ) 
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Therefore, by Cauchy-Schwartz and (47) 

Vu,p [ jn{u}r)(f^{a + br) 

J p ' 

This implies that the limit (56) is 0. 


dr 


< ^0, as p ^ 0. 


To take care of the other terms in /i, we choose 1 >> pi > p and consider 

^2 n / ^ 2 \ , p\2 


J*2 ri / 

„=lm=-n ^ “'Pi^ 


(a + bry 


dr. 


n—l m——n 

By (57), it is easy to see 

N 2 n 


rv2 

_ 1 ^ ^ J 


(a + try 


dr = 0. 


n—l m——n 


To see that the terms integrated on (p, pi) converge to the right hand side of (55), 
we first apply integration by parts 

rfn:p r {a+ br)^^^^ dr =<,(p-(l)^n(p) 

dp ^ 


rpi 


-d: 


n,p 


<p™'(a + br)bAn{r) dr 


where 


24„(r) := f h^\u!s){a + bs)'^s ^ ds. 

J r 

For 0 < n < -/V 2 , and r < pi << 1 


20r 


I \ "■+! rpi 


An{r)= - — ^r(n + 1/2) ( - ) / (a + 6s)^s (l + 0(s)) ds 


2^/^^ 

Therefore, as /? —;► 0 


UJ J 

r(., + i/ 2 )(/) 


n+1 2 

a 


n+1 


r-("+')(l + 0,^o(r"”)). 


Mp) = + Op-o(p)) 

2i/7r(n + 1) \UJ/ 


and by (57) 


(p™ (o + br)bAn{r) dr 


< 


r2 X 1/2 


npi u 

\Mr)\^ 


a + br 


dr 


^<p P^^^ —t 0 as p —?► 0. 


Together we have 
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From (57) and (44), we further have as p —>■ 0, 
(58) 

2zi/7r Jn{kuj)hn\ku:) - 'Hn{ku;)jn{ku}) 


d™ = 


implying 


( fjJ \ TL -\- 1 . 

2 ) (1 + Op^o(p)), 


r(n -I- 1 / 2 )^q ^^"^knjnikuj) 


p->-0 

Therefore, (55) is proved. 


1/2 

Mo 


Jn{kuj)hn\ku) - 'Hn{kuj)jn{kuj) 


kn{n + l)jn{kuj) 


Qn 


□ 


Before considering the tail term l 2 {p), let us define 

Bn(r) := f hl^\ujs){a + bs)'^s~^ ds, r > p. 
J r 

for n > N 2 (so we can use (78)), similar to A„(r), we have 


(59) 
and 

(60) 


Bn{p) — 


ir(n + 1/2) (2 


2y/n{n + 1) \a; 


n+1 


B„(r) = -^r(n + l/2) - 


2^/^T 


n+1 


p-("+')(l + Op^o(p)), 

p-("+i)(l + 0,^o(r-")) 


n -|- 1 


for n > N 2 . 

Lemma 4.4. 

(61) 


lim/ 2 (p)= E 

n=iV2 + l m——n 


1/2 

Mo 


Jn{kuj)hn\kuj) - 'Hri{kuj)jn{kuj) 


kn{n + l)jn{kuj) 




Remark 4.5. Along with Lemma 4-. 3, this shows that the limit of the normal 
component of the exterior field is some function (or distribution) times 5{r — 1). 
The smoothness of this function (the strength of the delta singularity) is estimated 
by the growth of the coefficient with respect to n » 1. For n > N 2 , by (77), (80) 
and (52), we have 


Si 


1/2 

,Mo 


Jn{kuj)hi)\kuj) - 'Hn{kuj)jn{kLo) 


kn{n + l)j„(fca;) 
-1 


1/2 { 2 n+l)T{n+l/ 2 ) ,, . 

_ Mo 2r(7i+3/2) 

M^+l) 2r(n?i/2) {^) 
for any M > 0. 


:SlK 


1 +o™(-)) <M 


kn{n + 1) ^ 


Proof of Lemma Taking into account the n-dependence, we have from (43), 
(38) and (44), 

(62) |<,|< r(n + l/2)r(n + 3/2)a2-"|g;r|(£)”^\ |<,| < |CI 

where the general constants associated to < are independent of n. 
























18 


3-DIMENSIONAL INVISIBILITY CLOAKING 


First to show 
(63) 


IL n ^ 

= E E ■■= / <p(0 dr, 

_AT- , 1 ^ Jo 


n—N2-\-l m——n 


by (62) and (77), we have for r G [p, 2] and n > N 2 


/ 2 \ 

<^/i(i)(a;r) +<^j„(a;r) < r(n+ 1/2) (^—j |C|p”+V-("+i). 

By the estimates (52) for and < N 2 , 


\hll\kujri)\ > r(n-h 1/2) 


kujr\ 


n +1 


we have 


Kp(r)| dr < ,52r(n + 1/2) (^—j |C|p”+' J K(a + &r)|r 

<M (l + n)-'^r/+V'/', 

implying that the series is uniformly convergent. So (63) is valid. 

To show 


-”-2 dr 


(64) 

lim 

p—^■O 

E E 


= E 

E lim^np 

o-yO 


n—N 2 -\-l m— — n 


n—N 2 -\-l m— — n 

we write for 

n > N 2 






'^Zp=Slr,zJ Mu^r)p] 
J p 

/(a + 6r)(a 

+ br)'^r ^ dr 



+ Sld^,p f h^n'>{ujr)p^{a + br){a + brfr ^ 


- 

^ l,n,p ' ^ 2,n,/ 

9' 



By (62) and 

(77), we 

have 




1 ^ l,n,pl 

< ^nlC 

|r(n + l/2)( 

p ^ 
2kuj ^ 

n+1 

1 / -”■ 
J P 

'Vn (a + &r)IO 


< 

r(n + l/2) / 
' 2n - 1 V 

' P 
~2ku! 

\n+l 

) (2" + 

0(p"-i/2)) 


< 

r(n + l/2) / 
n V 

. kuj/ 

n+1 



where the general constants are independent of n > N 2 and p. The right hand side 
is summable with respect to n uniformly in p, using (52). Moreover, it converges 
to 0 as p —T' 0. Therefore, 


(65) 


p™ E E ^F".p E E 

n=Ar2 + lm=— n n=N2 + l^= — 


For the second term, by the integration by parts, we have 

= Sld^,p |^“(l)i?„(p) - I' <'(a + br)bB^{r) drj . 
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Combining (62) and (77), we obtain 

. 2 r(n-b 1/2) 1 




< 


ri+1 \h^^\ku})\ 
< i^2^;r/i(i)(fcoi)<^:/(i)i 


152^/1(1) (fcu;)^™(l) I 


which is summable by that of (50). Also, by (62), (60) and (51), one has 


q2jm 

^n^n,p 


(/7™ (a + br)bBn{r) dr 

n +1 


< s'- 




PN "+1 r(n-b 1 / 2 ) 


n -b 1 


n +1 


p 


-(n+l/ 2 ) 


s s;ici'-- "// -■ (/j) (i + n)-<««)y/> 

By Lebesgue dominated convergence theorem, we have 


J“„ E E E E E E 

n—N2 + l'f^— — 'n n—N2-\-l'^——'ri n=A/’2 + l 7n=—n 

where 

:= hni <+„(p) 

p->0 


2ii/7r Jn{kuj)hn\ku:)-'Hn{kuj)jn(kuj) 


( /.I \ n +1 


r(n + 1/2)^q ^^'^knjn{kuj) 

by (58), which proves the lemma by (59). 


□ 


Remark 4.6. Using the above analysis, we can calculate the limit of the tangential 
components of the fields when we are cloaking an active source. Thus, one can 
obtain explicit representations for the the limits of right hand sides of equations 
(36) and (37). For example, we have 

T(]l{k,uj, J) := Wrn^Mkiv) + qffHn^ku:) = filf^Mkuj) + g™H„(fccc) ^ 0 

p—^■O 

where /3)/g = (53). Notice the convergence is uniform in n. 

Therefore, we have 

oo n 

( 66 ) lim f X EfloB, = +^ ^ SUTi]lVf^ + 

n—1 m——n 

where 

Ti%{k,u},J) := lini a™ j„(fcw) -bp™+^)(fcw) 

’ p—>-0 

Similarly, we can obtain the tangential magnetic boundary condition at the interface 
as well. 
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Appendix A. Spherical harmonics and Bessel functions 


Our arguments rely heavily on expanding the EM fields into series of spherical 
wave functions. To that end, we introduce for n G and to S Z, 

:= V X {xjn{uj\x\)Y^{x)}, 

Nn,Jx) := V X {xhl^\uj\x\)Y;^{x)}, 

where w S K and x = xl\x\ for x G Here, Y^{x) are spherical harmonics and 
hn\t) ;= jn{i) + Wnit) with jn(t) and yn(t), for t G M, being the spherical Bessel 
functions of the first and second kind, respectively. The following facts about these 
functions are useful in our estimates. 

Set Sn = \/n{n + 1). Define 

(68) Jn{t) := jnit) + Unit) := h^n\t) + 

where and hn'^ are the derivatives of j„ and hn'^ . We introduce the vector 
spherical harmonics 

1 


(69) [/^(x) := ^Grad ^^(x), ^^(x) := u x 

where Grad denotes the surface gradient. They satisfy 

(70) i X yr(i) = -G™(i), i X C/r(i) = Hr(f). 

Moreover, we can rewrite and as 

m:;^jx)= - s„j„{co\x\)v:^ixi 

-5„/iW(a;|x|)Hr(i)- 

Moreover, we have 

V X M™^(x) =Sn\x\-^MLu\x\)U:;^{x) + Sl\x\-^Uu:\x\)Y:rim, 

V X Af™^(a:) =S'„|a;r^'H„(w|ai|)G™(x) + (a;|a;|)y„'"(i)i. 

Alternatively, we also have by (69) 

V X M-Jcr) =1x1-1 J'„(c.|ai|)Vy7(x) + 52|a:|-ij„(a;|x|)yr(i)i, 

V X A^™^(x) =|xri'H„(a;|x|)Vy'™(x)+ S'^|xri/i(,i)(a;|x|)y'“(x)x. 

It is easy to see that 

x X M“^(x) = 5„j„(u;|x|)17;"(x), x x Ar™^(x) = 5„/iii)(c.|x|)C/™(x) 
X ■ M ^^^{ x ) = X ■ Af™^(x) = 0, 

and 

X X (V X M-^(X)) = Sn\x\-^Jn{uj\x\)V;^ix), 

X X (V X Ar™^(x)) = S'„|x|-i'H„(w|x|)E™(x), 


(71) 


(72) 


(73) 


(74) 


(75) 

(76) 


X • (V X M™^(x)) = Sl\x\ ij„(w|x|)y™(x), 
X- (V X iV-Jx)) = A2|x|-i/iW(a;|x|)yr(x). 


The spherical Bessel functions are given by 

jn{t) = V7r7(^J„+i/2(t), yn{t) = \/7r/(2t)y„+i/2 (f) 
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where Jn+i/ 2 {t) and F„+i/ 2 (i) are the standard Bessel functions. More specifically, 


, , sint . (ii, N sint cost 


Let r(n -I- 1/2) := \/tt. From their series representations, we obtain that 

(77) 


Jn(t) 


t 

2r(n-h3/2) V2 

in the following sense 




n+1 


for n » t, 


(78) 


jnit) = 


{l + 0{l/n)) 


2T{n + 3/2) \2^ 
as n —>■ oo, uniformly for t on a compact subset o/K, 

n+1 


'■•‘"W =(l) (1 + 0(1/")) 

as n —>■ oo, uniformly for t on a compact subset of (0, oo). 


and for each n > 0, 


(79) 


jn{t) = 


2r(n-f3/2) V2 


Zy TT 


(l + Ot^o{t)') , 


9 \ n+1 

-j (l + Oi-s.o(i)), 


where \Ot^o{t)\ < C{n)t as t —)■ 0 for some constant C{n) >0 depending on n. Due 
to (78), one can easily obtain a uniform C independent of n, hence replace On{t) 
by 0{t). (Notice that such uniformity for t small is corresponding graphically to 
the spreading out shape of jn(t) and yn{t) with respect to n, i.e., less oscillatory 
for larger n.) In the same sense, we have when n » t, 


(80) 


Jn{t) 


y/TT{n -I- 1) 

2r(n-h3/2) yj ’ 


H„(t) 


.r(n-f l/ 2 )n /2 

* [tj 
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